The Bimodal Ising Spin Glass in dimension three : Corrections to scaling 
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Equilibrium numerical data on the three dimensional bimodal (±J) interaction Ising spin glass 
up to size L — 48 show that corrections to scaling, which are known to be strong, behave in a 
non-monotonic manner with size. Extrapolation to the infinite size thermodynamic limit is difficult; 
however the large L data indicate that the ordering temperature T c lies significantly higher than the 
values which have been estimated from previous numerical work limited to smaller sizes. In view of 
the present results it is at the least premature to conclude that the three dimensional bimodal and 
Gaussian Ising spin glasses lie in the same universality class. 
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I. INTRODUCTION 

Renormalization group theory (RGT) provides an ex- 
planation of the physical origin of the critical exponents 
and of the universality classes for standard second order 
transitions which is one of the most remarkable achieve- 
ments of statistical physics. The universality rules state 
that the critical exponents depend only on a small num- 
ber of basic parameters^, essentially the dimension of 
space D, the range of interaction, and the number of 
order parameter components n. Physically, the critical 
parameters should not depend on the details of the short 
range interactions. The few known and well understood 
exceptions to universality concern mainly rare marginal 
cases, such as certain regularly frustrated spin systems 
in two dimensions where critical exponents vary contin- 
uously with the value of a control parameter. It should 
also be kept in mind that if the specific heat exponent a is 
positive for pure ferromagnets, disorder induces another 
universality class^. 

The glass transition in structural glasses has long 
raised fundamental questions. It has been hoped that 
the study of spin glass (SG) transitions would give insight 
into the basic physics of glassy ordering. In particular, a 
natural assumption to make is that universality rules hold 
in SGs, with exponents which are different from those of 
ferromagnets but which do not depend on the detailed 
form of the interactions between spins. 

For technical and conceptual reasons the SG transi- 
tion should be much easier to study than structural glass 
freezing but it has turned out that SG transitions also 
pose major problems. It has proved impossible so far to 
extend RGT satisfactorily to dimensions below the upper 
critical dimension, so information on transitions comes 
from experiment and the large scale numerical simula- 
tions which can be carried out on SGs. Numerous numer- 
ical studies on Ising Spin Glasses (ISGs) have concluded 
that universality holds also for these systems. However, 
determining transition temperatures and critical expo- 
nents to high precision at SG transitions through simula- 
tions remains a notoriously difficult task. Problems that 



are intrinsic to the numerical simulations include ago- 
nizingly slow equilibration rates near criticality and the 
need to average over a large number of microscopically 
inequivalent samples. It is helpful to analyze the numer- 
ical data using appropriate scaling variables and scaling 
rules^, and above all it is essential to allow adequately 
for corrections to finite size scaling (FSS). 



Here we address once again the case of the Ising 
spin glass with bimodal (±J) interactions in dimension 
three (3D), a system which has already been the sub- 
ject of many numerical studies. Recent careful measure- 
ments^^, where samples of sizes up to L — 24, 20, 24 
and L = 28 respectively were studied, all give estimates 
of the critical temperature consistent with T c ~ 1.12(1). 
In the present work equilibrium measurements have been 
made for sample sizes up to twice as large as those stud- 
ied in most of the previous work. The data demonstrate 
that in this system the influence of finite-size corrections 
to scaling extends to very large sizes, strongly affecting 
the thermodynamic limit estimate of the ordering tem- 
perature and therefore those of the critical exponents in 
this system. An analysis of these equilibrium data allow- 
ing for the corrections suggests a higher critical tempera- 
ture, compatible with estimations from a "dynamic" ap- 
proac h 8 i 9 i 10 . As a corollary the critical exponents will be 
significantly different from those obtained with T c ~ 1.12, 
as in Ref. 4,&iQ. 



Our principal conclusion is that the present large size 
data for the 3D Bimodal ISG taken together with pub- 
lished data for the 3D Gaussian ISG where corrections 
appear to be much weaker (Ref. @) cannot be read as ev- 
idence that the two systems lie in the same universality 
class. 
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II. MEASUREMENTS AND ANALYSIS 

We measure as usual the reduced spin-glass suscepti- 
bility x(T, L) defined as 



(1) 



where (• • ■) represents thermal and disorder average, and 
the upper suffix denotes two real replicas of the model 
with identical bond configurations. Using the spin over- 
lap given by 
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the spin-glass susceptibility is also expressed as \ = 
L 3 (q 2 ). We also measure finite-size correlation length^, 
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where x(k) is the Fourier transformed spin- glass sus- 
ceptibility with wave vector k and k m j n = (2vr/i, 0, 0). 
This becomes identical to the second moment correlation 
length in the thermodynamic limit L ^> £(oo, T). 

Published estimates of critical exponents on specific 
ISGs have varied widely (see the list in Ref. [(J) and it is 
of interest to examine the reasons for this dispersion. In 
most work FSS analyses are carried out on measurements 
made at equilibrium on small to moderate sized samples; 
the sizes L, the number of independent samples, and the 
equilibration times were necessarily restricted by com- 
puter time limitations. The ordering temperature T c is 
generally estimated by the intersection temperature of 
curves for the Binder parameter given by 
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or of the normalized correlation length £(L, T)/L as these 
two observables become independent of size at the order- 
ing temperature, to within FSS corrections. This caveat 
is vital because for the sizes which can be equilibrated 
in practice FSS corrections arc often not negligible and 
it is hard to extrapolate precisely from data on a set of 
finite L samples to the thermodynamic (infinite L) limit 
because, in contrast to the situation in standard ferro- 
magnets, there are no firm guidelines as to the form of 
the corrections at and close to criticality. We will return 
to this crucial point below. 

Once T c has been estimated the critical exponents v 
and rj (with 7 = (2 — 77)1/) have generally been esti- 
mated through a standard FSS analysis on the Binder 
parameter, the normalized correlation length, and the 
reduced susceptibility x(L,T) : g{L,T) = F g [L x / v (T - 
T c )], £(L,T)/L = F(\L 1 / V (T — T c )] and X (L,T) = 
L 2 -^F x [L l / v {T - T c )]. These rules, which have been 



taken over directly from the usage in ferromagnet stud- 
ies, implicitly assume that the appropriate temperature 
scaling variable is [(T—T c ) /T c ] . In a ferromagnet the cou- 
pling strength is J and so reduced field and temperature 
scales are H/J and T/J. In a SG the effective coupling 
strength parameter is (J 2 ) not J and it is the non-linear 
susceptibility, not the linear susceptibility, which diverges 
at criticality. Hence the appropriate "field" and "temper- 
ature" variables are H 2 / (J 2 ) a nd T 2 /(J 2 ) rather than 
H/J and T/J (see Refs. [l3IlMi3Lllb ■ We have shown 3 
that it is useful for SGs to employ "extended scaling" 
rules which (writing 0=1 /T) make the leading approx- 
imations 



£(oo,/3)~/3[l-(/3//3 c ) 2 



and 



x(oo,0)~[l-(/3//3 c ) 
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These reduce to the standard forms in the limit (T — 
T c )/T c <C 1, and they also lead to the functionally cor- 
rect forms in the high temperature limit. They become 
exact at all T > T c in the high dimension limit, but 
correction terms are to be expected in finite dimensions. 
When these scaling expressions are used, estimates for 
the exponents derived from scaling data on the different 
observables over wide temperature ranges are consider- 
ably more consistent than when the traditional rules are 
used to analyze the same data sets (see Ref. where the 
two analyses are compared). The traditional expressions 
reduce to the SG expressions very close to criticality, but 
if they are used to scale data covering a wider range of T 
they tend to strongly bias estimates of v. Consequently, 
previous works did not provide a consistent estimation 
of v that should be independent of physical quantities. 
This remark in itself explains why estimates of v in older 
publications were systematically very low compared to 
more recent values^. 

There remains the basic technical problem of deter- 
mining of an adequate sample size such that corrections 
to FSS have become negligible. Sizes L cannot be in- 
creased at will because equilibration at and near T c be- 
comes rapidly more difficult with increasing L, leading 
to escalating cost in computer time. We have chosen a 
strategy which sidesteps this problem by appealing to the 
finite-size ratio scaling approach introduced by Carac- 
ciolo et al — , and already applied to the 3D bimodal 
ISG by Palassini and Caraccioloii. For each observable 
Q(L, T) the ratio Q(sL, T)/Q(L, T) is plotted against the 
normalized correlation length £(L,T)/L. In practice we 
scale with s = 2 and so plot the ratios Q(2L, T)/Q(L, T), 
with Q = £, x an d <?• In the limit of sufficiently large L 
where finite size corrections to scaling have become negli- 
gible, for each observable Q(L, T) all the scaled data must 
fall on a single curve characteristic of the universality 
class of the system^. In this limit £(2L, T c ) /£(£, T c ) and 
g(2L,T c )/ g(L,T c ) are equal to 2 and 1 respectively when 
£(L,T)/L is equal to its critical value. Instead of only 
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using the intersections of £(L,T)/L and g(L,T) curves 
at and near T c to judge whether the FSS-correction-free 
limit has been reached, we can use the entire family of 
scaling curves. The important point is that if already the 
scaled data do not lie on a single i-independent curve at 
temperatures somewhat above T c , then finite size correc- 
tions to scaling are still important for the sizes studied 
and the thermodynamic limit has not yet been reached. 
Because of the high value of the dynamical exponent 
z{T c ) in ISGs (z(T c ) is typically ~ 6 in 3D) it is techni- 
cally much faster to equilibrate large samples down to say 
1.2T C so as to judge if corrections to scaling have become 
negligible rather than going to temperatures down to and 
below T c . The only drawback of this approach is that one 
needs to make an extrapolation in temperature (as well 
as in size) in order to estimate T c , the critical exponents, 
and the large size universal parameters !;(L,T C )/L and 

g(L,T c ). 

Cubic lattices of linear size L = 3 to L = 48 with 
bimodal random interactions with equal probability and 
conventional periodic boundary conditions were equili- 
brated using the exchange Monte-Carlo technique^. The 
number of independent samples studied were : 4000, 
4000, 4000, 8192, 8192, 6144, 1024, 2688 and 512 for 
L = 3, 4, 6, 8, 12, 16, 24, 32 and 48, respectively. Sam- 
ples up to L = 16 were equilibrated down to T = 0.9; in 
order to avoid the difficulties associated with excessive 
equilibration times at the lowest temperatures for large 
samples with L = 24,32, it was decided to equilibrate 
only down to T = 1.1 and 1.2 respectively. Two inde- 
pendent batches of 256 L = 48 samples were equilibrated 
down to T = 1.3 and to T — 1.4. As discussed above, for 
these large sizes this approach avoids the region below 
the ordering temperature where equilibration is particu- 
larly hard to achieve. The observables recorded were the 
spin-glass susceptibility x(L,T), the finite size correla- 
tion length £(L,T), the Binder parameter g(L,T), and 
the energy per spin e(L,T). Equilibration was checked 
by studying the dependence of the observables on the 
measurement time. Wherever direct comparisons could 
be made, i.e. for sizes up to L = 24, there was excel- 
lent point by point agreement between the present long 
time anneal values and equilibrated data from an inde- 
pendent study (Ref. [6| and Helmut Katzgraber, private 
communication), confirming that equilibrium had been 
reached. 

The data were analyzed using the Caracciolo et a££ fi- 
nite size ratio scaling approach. Figures [T] and [5] show the 
ratio scaling plots for the correlation length £(L,T). In 
order to be able to visualize the non-monotonic behavior 
we present the data in two scaling plots. It can be seen 
that for small L the scaling curves move systematically 
from left to right with increasing L; then as L is increased 
further beyond L ~ 24 the curves begin to move to the 
left again. By inspection the present data demonstrate 
that in the 3D bimodal ISC there are significant finite- 
size corrections to scaling up to and probably beyond the 
largest sizes which we have studied; the observed non- 
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FIG. 1: (Color online) Scaling plot for correlation length ra- 
tios £(2L,T)/£(L,T) against £{L,T)/L following Ref. [H. 
Top to bottom L = 3, 4, 6, 8, 12 cyan, blue, green, red, black. 
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FIG. 2: (Color online) Scaling plot for correlation length ra- 
tios £(2L,T)/£(L,T) against £(L,T)/L following Ref. [H. 
Top to bottom L = 24, 16, 12 black, red, green. 



monotonic behavior of the scaling plots as functions of L 
implies that there are at least two important and conflict- 
ing correction terms. The critical normalized correlation 
length (£(L,T)/L) C corresponds to the x coordinate of 
the intersection of the limiting large size curve with the 
horizontal line y = 2. Data taken only up to L ~ 24 give 
the impression that the large size limit has been reached 
by this size, and that (£(L,T)/L) C ~ 0.65 (see Katz- 
graber et al and Hasenbusch et al&£). However when 
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FIG. 3: (Color online) Scaling plot for reduced susceptibility 
ratios x(2L,T)/x(L,T) against £(L,T)/L following Ref. Q]| 
Top to bottom L = 3, 4, 6, 8, 12 : cyan, blue, green, red, black. 



FIG. 5: (Color online) Scaling plot for Binder parameter ra- 
tios g(2L,T)/g(L,T) against £(L,T)/L following Ref. [H. 
L — 24, 16, 12, 8,6,4 : black, red, green, blue, cyan, magenta. 




FIG. 4: (Color online) Scaling plot for reduced susceptibility 
ratios x(2L, T)/x(L, T) against £(L,T)/L following Ref. EE 
Top to bottom L = 24, 16, 12 : black, red, green. 



FIG. 6: (Color online) Scaling plot of the Binder parame- 
ter g{L,T) against £(L,T)/L. L = 48,32,24,16,12,8,6,4 : 
black, red, green, blue, cyan, magenta, yellow, dark yellow. 



data up to larger L including L = 48 are taken into ac- 
count, it is clear that one needs to go far beyond L ~ 24 
to reach the limit where FSS corrections have become 
negligible; the data indicate that the true large size limit 
for (£(L,T)/ L) c will lie below 0.65, and will correspond 
to an ordering temperature T c in the infinite size limit 
significantly higher than recent estimates^£i££ 

Figures [3] and Q] show the analogous ratio scaling plot 
for the SG susceptibility x(L,T). Again for small L the 



curves move systematically to the right with increasing 
L, and then begin to move to the left for larger L. Once 
more the position of the infinite L limit curve is diffi- 
cult to judge but it would correspond to a value for the 
exponent rj less negative than that suggested by recent 
estimates^. Figure [5] shows a ratio scaling plot for the 
Binder parameter g(L, T). Again similar behavior can be 
observed. It is of interest to note that when the same data 
are represented in the form of a plot of g(L,T) against 
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FIG. 7: Correlation length ratios £(2L,T)/£(L,T) against L 
for £(L,T)/L = 0.393. 



FIG. 8: Susceptibility ratios x(2L,T)/x(L,T) against L for 
£{L,T)/L = 0.393. 



£(L,T)/L as in Fig. [6l for L > 4 all the points lie on a 
single curve within the present statistics. This form of 
scaling plot is very insensitive to corrections to scaling 
because g(L,T) and £{L,T)/L are closely correlated. 

Comparisons have been made between different spin 
glasses using this form of scaling plot (e.g. Ref. H) 
because systems lying in different universality classes 
should not have the same scaling curves. However simple 
inspection of figures [U [21 [3] and 0] shows that g(L,T) re- 
mains an almost unique function of £(£, T)/L even when 
strong finite size corrections to scaling are obvious in the 
ratio scaling plots for each of the two parameters taken 
individually. In other words the corrections to scaling for 
g(L,T) and for £(L,T)/L are strongly correlated. Be- 
cause this particular form of scaling is very insensitive to 
corrections, using such plots to judge if two systems lie 
in the same universality class or not would require data 
of extremely high statistical precision. 

As an illustration, ratio data are shown as functions 
of L in figures [3 [8] and [9] for one fixed normalized 
correlation length value, £(L,T)/L = 0.393. The fig- 
ures demonstrate the non-monotonic behavior, and also 
the difficulty in extrapolating to infinite L when aim- 
ing to estimate the limiting ratio values, even at tem- 
peratures well above criticality. If nevertheless we as- 
sume that the £(48, T)/£(24, T) ratio curve is close to 
the infinite size limit, and we extrapolate the curve to 
£(48,T)/£(24,T) = 2 we obtain T c ~ 1.175, which in 
view of the extrapolations we have been obliged to make 
can be considered as a lower limit to the true value of T c . 
We recall that a quite independent "dynamic" estimate 9 
gives T c ~ 1.19. This value is consistent with the present 
equilibrium analysis. 

It is important to envisage potential sources of sys- 
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FIG. 9: Binder parameter ratios g(2L,T)/ g(L,T) against L 
for £{L,T)/L = 0.393. 



tematic error in the present data. Even though stringent 
precautions have been taken to assure full equilibration, 
as always with spin glass simulations a possible source of 
error which must be kept in mind would be incomplete 
equilibration. As stated above, the data for sizes up to 
L = 24 arc in very good agreement with those from reli- 
able independent studies (Ref. and Helmut Katzgraber, 
private communication). It should be noted that a lack of 
equilibration for the larger samples only [L = 48 and L = 
32) would reduce £(2L, T) and x(2L, T) for L = 24 and 
L = 16 while leaving £(L,T),x(L,T) and £(L,T)/L un- 



6 



altered, and so would artificially move the scaling curves 
of Figures 2 and 4 downwards. The observation that 
the curves for £(48, T)/£(24, T) and £(32, T)/£(l 6, T) in 
these two figures lie well above the others thus cannot be 
an artifact arising from lack of equilibration at L = 48 
and L = 32. There remain of course purely statistical 
sampling effects due to the fact that at large sizes the 
number of samples studied is necessarily restricted. For 
the sizes L = 48 and L — 32 the only data with which 
comparisons could be made are those from the pioneering 
work of Ref. Il7l : these are fragmentary and recorded on 
far fewer samples than in the present work, but within the 
statistical errors there is agreement, particularly for the 
two L = 48 points from 64 samples recorded by Ref. 0- 
For the 3D Gaussian ISG, corrections to scaling are 
much very weaker than in the bimodal case^, so it can 
plausibly be assumed that the estimates for T c and for 
the critical values of £(L,T)/L and g(L,T) from the 
largest L Gaussian ISG data- can be taken as being es- 
sentially equal to the infinite size thermodynamic limit 
values. Unless strong corrections appear when still larger 
L Gaussian samples arc studied, the Gaussian critical 
parameters in the thermodynamic limit are not equal to 
those for the thermodynamic limit bimodal ISG derived 
as outlined above, implying that the two systems do not 
lie in the same universality class. 

III. FINITE-SIZE CORRECTIONS TO SCALING 

Clearly it is essential to allow for corrections to scaling 
in order to obtain reliable estimates of the critical param- 
eters. At criticality FSS expressions including correction 
terms can be written as (see for instance Ref. [Til ) 

X (L, f3 c ) = C X L 2 -^{1 + ai L-" + a 2 L-"* + . . .] (7) 

where u> and u>2 are the first and second correction ex- 
ponents. For standard ferromagnets these exponents 
are rather accurately known from RGT; in dimension 3 
lu ~ 0.51 and u)% ~ 0.82 in pure systems and u) ~ 0.3 and 
uj2 ~ 0.8— for diluted ferromagnets. Unfortunately there 
are no equivalent RGT guide-lines for the values of the 
correction exponents (and a fortiori for the strengths of 
the correction terms) in ISGs. The non-monotonic varia- 
tion of the positions of the scaling curves with increasing 
L in the 3D bimodal ISG implies that there are two sig- 
nificant contributions (at least) with opposite signs to 
the finite-size-scaling corrections, coming from to and L02 
terms. 

For finite-size scaling with corrections near T c , Cal- 
brese et al— give the ansatz 

£(2L, T)/£(L, T) = F(L/£(L, T)) + L^G(L/^L, T)) 

(8) 



and the equivalent equation for \i where F and G are 
scaling functions and u> is the leading non-analytic cor- 
rection to scaling exponent. Unfortunately the present 
data are manifestly insufficient to attempt fits to anal- 
ogous equations with two correction terms, so we can 
make no sensible estimate of the correction exponents in 
the 3D bimodal ISG. 

It can be noted that in the 2d Bimodal ISG (which 
orders only at T = 0) limiting behavior at criticality is 
only attained above L ~ lOO^i. 



IV. CONCLUSION 

Equilibrium numerical data on the three dimensional 
bimodal Ising spin glass for samples up to size L — 48 
show that corrections to finite-size scaling are unexpect- 
edly strong and non-monotonic. Taking into account the 
influence of these corrections, it can be concluded that 
the large size limit behavior and the corresponding ther- 
modynamic limit critical temperature and critical param- 
eters are significantly different from the values derived 
from previous equilibrium simulation data where smaller 
maximum sizes were studied. Estimates of critical pa- 
rameters from the present data are only indicative due to 
the intrinsic problem of extrapolating reliably to infinite 
size, but they are consistent with those obtained from 
a "dynamic" method£&i£. The strong corrections to fi- 
nite size scaling must be fully taken into account before 
a valid comparison can be made between the thermody- 
namic limit behaviors for the bimodal and Gaussian 3D 
ISGs. The data as they exist cannot be considered to 
show that the two systems lie in the same universality 
class. 
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